function P=AVP_step(G,P,k)
tspan=[G.t,G.t+G.h];
y0=P.Vars;
[t,y]=ode23(P.ode_func,tspan,y0,[],G,P,k);
P.trace=y;
P.Vars=y(end,:);
%~ P=feval(P.calc_pos_ang_func,G,P,k);

